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I present a theory-guided review of RHIC data, arguing that they provide strong ev- 
idence for formation of a thermalized quark-gluon plasma at RHIC. Strong radial flow 
! reflects high thermal pressure in the reaction zone. Large elliptic flow proves that the 
pressure builds up quickly and the system thermalizes on a very short time scale of 
<lfm/c. The observed hadron yields are consistent with statistical hadron formation 
from a quark-gluon plasma, followed by immediate chemical decoupling due to strong 
radial expansion. The observed suppression of jets appears to confirm the predicted large 
energy loss suffered by hard partons moving through a quark-gluon plasma; more work 
is required to quantitatively understand this effect. Source size measurements using two- 
particle correlations do not seem to fit into this picture; the origin of this discrepancy 
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("HBT puzzle") is presently not understood. 
1. THE LITTLE BANG: QGP AT RHIC 



A quark-gluon plasma (QGP) is a thermalized system of deconfined quarks, antiquarks, 
and gluons. As such it has thermodynamic pressure P = (? s e where e is the energy density 
and c s is the sound velocity. Perturbative QCD gives c 2 s = l/3 at leading order, and 
lattice QCD confirms [ Q that for T > 2T C (where T c « 170 MeV is the critical temperature 
for color deconfinement in QCD [ p]]) about 80-85% of this value is reached. In heavy 
ion collisions the reaction zone is surrounded by vacuum with zero pressure; thus, if 
the collision fireball contains a QGP, the pressure gradient near the surface will lead to 
collective expansion ("flow"), in particular transverse to the beam in which direction the 
nuclear matter was initially at rest. Therefore, collective transverse flow is an unavoidable 
consequence of QGP formation, and the data must show it if a QGP has been created. 

The converse is not necessarily true unless the observed flow is so strong that only a 
QGP could have generated sufficient pressure over a sufficiently long time to create it. In 
order to assess whether the latter is the case one exploits three facts: 

(i) Due to incomplete stopping of the two nuclei at high collision energies, the reaction 
zone expands quickly in the beam direction, thereby rapidly cooling and diluting the 
matter inside; hence the pressure has only a limited time interval available to generate flow 
in the transverse directions. This time is the shorter the lower the initial energy density; 
one can therefore establish an upper limit for the amount of "radial" transverse flow that 
can be generated if the initial energy density (at the time of approximate thermalization) 
was never significantly above the critical value for color deconfinement (about 1 GeV/fm 3 ). 
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(ii) By energy conservation, the initial energy density is related to the final transver- 
sally emitted energy per unit rapidity, dE?/dy, via geometric factors and the time of 
thermalization [|2[]: the shorter the thermalization time, the higher the initial energy den- 
sity at fixed (measured) final dEx/dy. The Bjorken formula [§] £(r t h or m) = ^ — ^ay^ff 
is actually an underestimate of the initial energy density since it neglects longitudinal 
work done by the pressure during the expansion [ ||. The measured values for dE^jdy at 
RHIC [ |7j are consistent with subcritical initial energy densities only if one assumes that 
thermalization takes at least r t h er m > 5 fm/c. 

(iii) On the other hand, this thermalization time scale can be constrained by the mea- 
sured elliptic flow in noncentral collisions [ ||. In such collisions, the nuclear reaction 
zone is initially deformed in coordinate space: a cut transverse to the beam direction 
looks like an almond whose longer side points perpendicular to the reaction plane. For a 
given collision centrality, extracted e.g. from the measured total multiplicity dN/dy, this 
deformation can be calculated from the overlap geometry [ [5| . Thermalization then leads 
to pressure gradients which are anisotropic and larger in the short direction of the almond, 
causing a faster growth of the transverse flow into the reaction plane than perpendicular 
to it. This mechanism leads to an anisotropy in the final transverse momentum distribu- 
tion, making it flatter in the direction of the impact parameter vector than perpendicular 
to it [ ||. The elliptic flow coefficient v 2 , defined as the second harmonic coefficient in 
an azimuthal Fourier decomposition of the measured transverse momentum spectrum, 
quantifies this anisotropy. Without reinteractions among the produced particles, no such 
momentum anisotropy would arise; rescattering in the hot fireball matter is required to 
transfer the initial spatial anisotropy into momentum space. Microscopic calculations 
[ m show that, for a given initial spatial deformation, the generated amount of v 2 is a 
monotonic function of the mean free path (or the product of density and scattering cross 
section) in the fireball; the maximal value is reached in the hydrodynamic limit of zero 
mean free path [ |S| . If thermalization is delayed, the initial spatial deformation disappears 
spontaneously by transverse free-streaming [ [jlj, thereby reducing the maximal possible 
momentum anisotropy. Thus, the measured elliptic flow at given collision centrality pro- 
vides immediately an upper limit for the thermalization time scale, with larger i> 2 values 
corresponding to smaller values for r t h C rm- In the hydrodynamic limit the extracted r t h e rm 
depends parametrically on the stiffness of the equation of state (velocity of sound c 2 ), 
with smaller sound velocities resulting in shorter thermalization times at fixed v 2 [ [TOf] . 
The weakest upper limit for r therm is thus extracted from the v 2 data by taking the largest 
reasonable value for the speed of sound, c 2 — |, corresponding to an ideal gas of massles 
quarks and gluons. 

Experimentally, the strength of the radial flow can be extracted from the transverse 



mass (m_i_ = Jm 2 + p±) spectra of a variety of different mass hadrons by fitting them 



to a common flow spectrum [ nj. At high m_i_3>m, the transverse flow flattens the 
thermal spectra by a simple common blueshift factor, and all spectra approach the same 
asymptotic slope from which the decoupling temperature and flow velocity cannot be 
separated. At low m± < 2m, however, the flow induces an even stronger flattening which 
increases with the rest mass of the particles. Heavier hadrons thus develop a visible 
shoulder at small m±, and this feature can be used to determine the average transverse 
flow velocity uniquely. An excellent systematic study of the shapes of the m±_ spectra of 
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7T ,p,p, A, A, H,H, f2,Q, and d from Pb+Pb collisions at beam energies of 40, 80 

and 160 A GeV at the SPS was recently presented by M. van Leeuwen at the Quark Matter 



2002 conference [ O] , showing perfect consistency of all these spectra with a two-parameter 



flow fit [ |TTJ with kinetic freeze-out temperature Tf = 120—130 MeV and average transverse 
flow velocity (P±) =0.4 — 0.5 (Tf and (j3±) are anticorrelated in these fits). Even the fl 
and fl show a clear flow shoulder and are consistent with these low freeze-out temperature 
and large flow values. This contradicts earlier conclusions drawn from simple exponential 
fits to fl spectra measured at higher m± values by the WA97 Collaboration [ [13|] which 
did not show clear evidence for the flow shoulder and which led to the suggestion [ [14| 
that fl and fl are too heavy and too weakly coupled to the expanding pion fluid to pick 
up much of the transverse flow created during the late hadronic stages of the collision, 
allowing them to decouple earlier, i.e. at higher temperature and smaller flow. The new 
data [ [12| show that the fl and fl fully participate in the flow and essentially decouple 



together with the rest of the hadron fluid. They also demonstrate the fallacies connected 
with characterising the spectra by a single slope parameter instead of comparing their 
entire non-exponential shape with a proper flow parametrization. 

The hadron data at RHIC are already almost as rich as at the SPS, but not yet all 
published in final form. Available flow analyses of pion, kaon and (anti)proton spectra [ 
|i~5| , |i~6l , |17|| give a similar range of freeze-out temperatures and average flow velocities as 
at the SPS, but in a different combination (either T f or (/3j_) is about 5-10% higher than 
at the SPS). This results in about 5% flatter pion spectra while the effect on the proton 
slope at small p± is much larger and at least 25%. Again the flow fits describe all available 
spectra very well up to transverse momenta of about 2-3 GeV (which covers more than 
99% of all hadrons), and only at higher p± one begins to see evidence for the power-law 
tails expected from hard QCD processes. Also at RHIC, preliminary fl and fl spectra [ 
|18[1 are consistent with the fl fully participating in the hadronic flow. A comparison with 
hydrodynamic predictions published in [ [0| works well if the fl is assumed to decouple at 
T f ~ 135 - 140 MeV, i.e. only slightly before the pions freeze out at T f = 125 - 130 MeV [ 
PH , but not if one assumes decoupling already at hadronization, Tf = T c = 170 MeV, where 
the hydrodynamic model has not yet developed enough transverse flow. 

The extracted flow velocities at or above half the speed of light demonstrate that the 
collision fireball indeed undergoes a violent explosion - the "Little Bang" . It appears to 
be impossible to obtain such large flow velocities without assuming initial energy densities 
and pressures well above the critical value for deconfinement. This is already true at the 
SPS [ [H]] but, as I show next, much more convincingly so at RHIC. 



A decisive measurement is the observation at RHIC, first made by STAR [ |22| and 



then confirmed by PHENIX [ ||, |4f and PHOBOS [|5]J, that the elliptic flow v 2 in non- 
central collisions is large and, for transverse momenta below p± ~ 2 GeV, almost exhausts 
the upper limit provided by earlier hydrodynamic predictions [ [K], [1^, |26| (see Fig. [I]). 
The agreement between theory and data requires that the hydrodynamic evolution starts 
no later than about lfm/c after nuclear impact [ [10], pl| (the successful predictions in 



[10| , |26l use Ttherm = 0.6 fm/c). Since the spatial deformation responsible for the creation of 



flow anisotropies quickly decreases, v 2 develops and saturates early in the collision [ |, |10| , 
long before hadrons decouple. Hence, the conclusion that the measured large v 2 implies 
early thermalization is not changed [ Wh if the somewhat unrealistic sharp "Cooper- Frye 
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freeze-out" used in the hydro dynamic approach is replaced by a proper kinetic treatment 
of the freeze-out stage [ |2"8fl . The agreement of ^(rl) out to transverse momenta of 
about 2 GeV, including the predicted [ [H| splitting of v 2 with the hadron rest mass (see 
Fig. |l]), shows that the bulk of the matter (i.e. > 99% of the emitted hadrons) behaves 
hydrodynamically. (One should note, however, that the measured rapidity dependence of 
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Figure 1. Transverse mo- 
mentum dependence of v 2 
for identified particles, 7r~, 
K~, p (left) and vr+, K + , 
p (right), as measured by 
PHENIX [ H. The 
solid lines represent the 
hydrodynamic prediction [ 
|T9"| for a freeze out tem- 
perature Tf = 120 MeV for 
(from top to bottom) ir, K, 
and p, respectively. 
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t>2 [ |25| cannot be reproduced by existing hydrodynamic models [ p9].) Microscopic mod- 



els share this success only if they approach the hydrodynamic limit by invoking unusually 
strong rescattering among the fireball constituents [ |30], |3~1"|] , far above the level con- 
ventionally expected from perturbative QCD. It is important to realize that, given the 
observed final multiplicity and transverse energy per unit rapidity, the maximum energy 
density in the fireball at r t hcrm ~ 0.6fm/c is about 25GeV/fm 3 in central collisions. Even 
after averaging over the transverse plane this is still more than an order of magnitude 
above the critical value for deconfinement and corresponds to about twice the critical tem- 
perature! Unless a completely different mechanism for creating the elliptic flow can be 
found, the conclusion seems unavoidable that a thermalized QGP at e > 10 e c and T > 2T C 
is created at RHIC which lives for about 5-7fm/c before becoming sufficiently dilute to 
form hadrons. Perturbative mechanisms [ |32| seem unable to explain the phenomenolog- 
ically required very short thermalization time scale, pointing to strong non-perturbative 
dynamics in the QGP even at or above 2T C . 

2. STATISTICAL HADRONIZATION: MEASURING T c 

The quark-hadron phase transition is arguably the most strongly coupled regime of 
QCD. Soft hadronization happens through a multitude of different channels and is there- 
fore most efficiently described in a statistical approach, as realized by Hagedorn more than 



35 years ago [ [33j . The microscopic processes are only constrained by local conservation 
laws (valid inside causally connected volume regions AV^) for energy, baryon number and 
strangeness. Maximizing the entropy subject to these constraints results in local thermal 
and chemical "equilibrium" distributions for the hadrons [ |34| whose local temperature 
T chem and chemical potentials fis,^s arise from Lagrange multipliers and reflect the lo- 
cal energy, baryon and strangeness densities at hadronization. This "equilibrium" is not 
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achieved kinetically (by hadronic rescattering) , but statistically (by interference of many 
channels). It does not require the prehadronic state to be thermalized, although very 
strong deviations from thermal equilibrium (e.g. at high p±) may survive the hadroniza- 
tion process. The statistical approach is only expected to work for soft hadron production. 

We know from lattice QCD that a hadronic equilibrium distribution at e > e c is unsta- 
ble against deconfmement, so statistical hadronization can only proceed once the energy 
density has dropped to £had = £ c ~ 1 GeV/fm 3 [|35|. The resulting equilibrium distribution 
will thus be characterized by T c h em — T c . Hadronic cascade models have shown [ [36[ that, 
at small net baryon density, the total particle density below T c is too small and the ex- 
pansion of the collision fireball is too rapid to maintain chemical equilibrium by inelastic 
hadronic collisions. Hence, the hadron abundances decouple directly at hadronization, 
and the chemical freeze-out temperature extracted from a statistical model fit to the 
measured abundances reflects directly the (de) confinement temperature, T c ^ em — T c [ [37 |. 
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Figure 2. Particle ratios 
from 130AGeV Au+Au 
collisions, measured by 
the STAR Collabora- 
tion at RHIC [ [18| and 
fitted by M. Kaneta 
to a thermal statistical 
model 
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Figure |2] [ [Tj| shows that the RHIC data confirm this expectation, just as it was previ- 
ously confirmed at the SPS [0: the chemical decoupling temperature extracted from the 
particle ratios is consistent with the critical temperature for deconfmement, T c ~ 170 MeV, 
extracted from lattice QCD simulations [ |IJ . Fitting the maximum entropy parameters to 
some of the more abundant particle yields one is able to reproduce all measured particle 
ratios, including those involving the rare multistrange (anti)baryons; this is quite impres- 
sive. Thermal freeze- out of the momentum distributions is delayed by strong quasi-elastic 
scattering via hadron resonances, such as7r + A^^A^7r + A^, tt + N^A^tt + N, 
tt + K — > K* — > 7i + K, etc., which do not modify the total measured yields of pions, 
kaons, and nucleons but continue to adjust their momentum distributions to the falling 
temperature and growing collective transverse flow until also resonance scattering ceases. 
The momentum spectra thus reflect a lower freeze-out temperature T f < T chem , as seen in 
section |l|. Note that only the total abundances (after all decays of unstable resonances 
are taken into account, e.g. iV^ ot = N n + 2N p + ]Va + N K * + . . .) freeze out at T chem while 
the fraction of pions stored in resonances still changes: due to their strong coupling to 
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the cooling pion fluid, the resonance abundances keep readjusting to the decreasing tem- 
perature, and the larger their pionic decay width the later they decouple. This implies 
that the measured K* /K, A/N, A* /A etc. ratios should generally be smaller than their 
chemical equilibrium values at T c h em = T c and, when themselves translated into a decou- 



pling temperature, reflect more closely the spectral temperature Tf than T c [ ^TJ. First 
results on resonance/ground-state ratios were reported at Quark Matter 2002 [|42j], and 
this prediction will be checked soon. 

3. JET QUENCHING: STRONG PARTON ENERGY LOSS IN THE QGP 

In 1982 Bjorken [ [E| suggested that fast partons travelling through a QGP might 
lose large amounts of energy by elastic scattering with the plasma constituents, resulting 
in the suppression of jets from the interior of the collision fireball in relativistic heavy 
ion collisions. Although the energy loss mechanism envisaged by him turned out to 
be ineffective, his qualitative prediction appears to be impressively confirmed by recent 
RHIC data. A good review of the theory of parton energy loss was recently given by 
Baier [ |4^|, and I refer to his talk for references. In a series of papers in the early 1990's, 



Gyulassy, Pliimer and X.N. Wang identified radiative energy loss as the dominant jet 
suppression mechanism, and in 1995-1997 Baier, Dokshitzer, Mueller, Peigne and Schiff 
(BDMPS) realized that an important effect on the energy loss rate results from multiple 
color interactions of the radiated gluon with the colored plasma constituents [ f45|. In the 



limit of an optically thick quark-gluon plasma they found [ [|5[ that the energy loss AE of 



the fast parton increases quadratically with the distance L travelled before escaping and 
hadronizing, 

AE«f£L>, with ^/<^i;|, (1) 

where a s is the strong coupling constant, p is the plasma density, \j? is the Debye screening 
mass for color electric fields in the plasma and A is the gluon mean free path. This was 
subsequently improved upon for plasmas with finite opacity by Gyulassy, Levai, and Vitev 
[ 16| and by Wiedemann [ |47]] , using an opacity expansion. For optically thin plasmas 



this leads to an important dependence of the energy loss AE on the energy E of the fast 



parton [ [46[ and to deviations from the quadratic path length dependence for small values 
of L. Both corrections reduce the predicted energy loss relative to the BDPMS result (|I|), 
in particular at "low" (SPS and RHIC) energies, but it remains still significantly larger 



in a QGP than in subcritical normal hadronic matter [ [BTj . With these improved results 
one may nurture the hope that parton energy loss and jet quenching, if observed, can be 
used as a quantitative probe of the density and its early time evolution in quark-gluon 
plasmas created in heavy-ion collisions. 

Two important observations at RHIC have moved this hope to the brink of reality. The 
first is the discovery of a suppression of high-p^ particle production in central Au+Au 
collisions at RHIC. When studied as a function of collision centrality, the production rates 
of charged hadrons at high p± 3> 2 GeV were found [ f|5| |5IJ not to follow the expected 
scaling with the number of binary collisions which characterizes hard QCD processes, 
but instead to continue to scale with A" part (the number of wounded nucleons), just as 
in the soft \ow-p± regime where phase cherence and the Landau-Pomeranchuk-Migdal 
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effect are known to suppress particle production. (There seems to be an intermediate p± 
region, 1.5GeV< p± <3.5GeV, where particle production grows a bit faster than N pait 
but more slowly than the number of binary collisions.) From the Glauber model the 
number of binary iViV collisions is known to scale with N p ^ t while iVp art is proportional 
to the volume of nuclear matter participating in the collision. The observed suppression 

1 /3 

or "quenching" of hard particle production by a factor iV^t is thus proportional to the 
linear size ("radius") of the nuclear overlap volume. It is qualitatively consistent with 
expectations from pQCD-based jet quenching calculations [ and inconsistent with 
pQCD calculations without jet quenching. It is seen for charged and neutral pions [ |50| 
but apparently not for high-p^ protons [ ^2] (at least not in the presently accessible p± 
range for proton PID); to what extent this "lack of hard proton suppression" is simply 
a reflection of the observed increase of the p/ir ratio with increasing pj_ [ Q which, at 
least up to pj_ ~ 2 GeV, is naturally explained by the stronger transverse flow effects on 
protons than on pions, is still unclear. 
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Figure 3. Ratio of Au+Au 
and p+p of correlated pairs 
(see [ Q for details) for 
small (squares) and large 
(back-to-back, circles) az- 
imuthal opening angles, for 
two different trigger p? in- 
tervals (left and right pan- 
els). The dominant system- 
atic error arises from the un- 
certainty in the elliptic flow 
v 2 at these pr- 

The second important observation, made by STAR [ |54|] and (not yet as convincingly) 
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by PHENIX 



], is that high-px particle production is indeed due to jets, and that the 
observed overall suppression of high-p^ particle production in central Au+Au collisions 
goes hand-in-hand with the complete suppression of the far-side partners of the observed 
jets. Perturbatively, hadron jets result from hard parton collisions leading to parton pairs 
with large, approximately balancing transverse momenta which fragment into pairs of 
back-to-back jets. In an angular distribution around the beam axis, jets appear in the 
two-particle correlation function as two peaks separated by 180°. By triggering on a high- 
p± particle and correlating it with other particles of p± > 2 GeV (in order to suppress 
the uncorrelated soft background), these pairs of peaks were indeed identified in pp and 
peripheral Au+Au collisions at 200 A GeV at RHIC [ [54j] , but in central Au+Au collisions 
only the correlation peak at small relative angles is found. Already in peripheral Au+Au 
collisions the far-side peak at 180° is much smaller than expected from a superposition of 
55j, and in central collisions it completely disappears after accounting 

151 



pp collisions [ p4 



(see Fig. 



for the collective azimuthal correlations generated by elliptic flow 

While attempts to quantitatively understand these observations have only just begun [ 
|5fJ , they suggest a very simple and intriguing picture similar to the one already proposed 
in 1982 by Bjorken [ 03]: fast partons formed inside the hot and dense collision zone 



s 



suffer severe energy loss and become part of the approximately thermalized, collectively 
expanding "soup" which emits particles at low p± < 2 GeV. Only if the primary hard 
collision happens near the surface of the reaction zone, the "outward moving" parton 
fragments to become a jet, very much like in pp collisions (upper set of symbols in Fig. |||), 
whereas the "inward moving" parton is more and more efficiently absorbed by the medium 
as the nuclei collide more centrally and the reaction fireball volume increases (lower set of 
symbols in Fig. ^). For central Au+Au collisions and p± < 8 — 10 GeV [ [54|, the "inward 
moving" parton doesn't make it to the other side with sufficient energy left to fragment 
into a recognizable jet. Jets are therefore only emitted from a relatively thin surface layer 
while the bulk of the fireball volume is opaque to jets. 

This very qualitative picture may appear naive but it nicely explains as a surface/volume 

1/3 

effect the missing factor N part in the observed scaling of high-p^ particle production. It 
also allows for a simple geometric estimate of the elliptic flow t> 2 at high p^, resulting 
from the anisotropic emission of surface jets from the spatially deformed overlap region [ 



57fl . As pointed out by Voloshin [ j58[ , this simple geometric estimate [ |57| appears to be 
consistent with the STAR data for v 2 at high p±. Of course, at very high p±_ one eventually 
expects the "inward moving" parton to emerge on the other side with sufficient energy 
left to form a jet. When this happens, its energy loss can be studied as a function of its 
path length through the medium, in order to check the energy loss formula ([I]) and its 
various published variants. Selecting events of fixed centrality (impact parameter), the 
path length can be controlled geometrically by the azimuthal angle between the jet and 



the reaction plane, and a new era of "jet tomography" [ 59] will begin 



4. THE RHIC HBT PUZZLE 

Unfortunately, space-time limitations do not allow for a detailed discussion of one piece 
of the puzzle which so far fails to seamlessly fit into the picture painted above: the 
two-particle Hanbury Brown - Twiss (HBT) correlation measurements. Since I reviewed 
these elsewhere (see [ |2D[ and references therein), let me be brief here. HBT correlations 
can be used to probe the size, shape and dynamical state of the source at hadronic 
decoupling. STAR and PHENIX have produced mutually consistent results showing that 
the HBT radii in the sideward and outward directions (i.e. perpendicular and parallel to 
the transverse emission vector) depend strongly on the transverse pair momentum and 
are almost equal to each other. Otherwise successful dynamical models such as those 
discussed in Sec. [I] fail to reproduce this strong transverse momentum dependence and 
consistently give smaller sideward than outward radii, in contradiction with the data. 
They also overpredict the longitudinal radius. The longitudinal and outward radii can 
be decreased by reducing the total lifetime of the fireball, but this would require even 
faster thermalization to accomodate the observed transverse flow, and the sideward radius 
remains too small. Equal outward and sideward radii not only require a very opaque 
source (such as the hydrodynamic ones with unrealistic sharp Cooper- Frye freeze-out), 
but also a short lifetime and a breaking of longitudinal boost invariance, at least in the 
decoupling process. While these problems are generally believed to reflect insufficiencies 
in the description of the late decoupling stage which will not affect our understanding of 
the early collision (QGP) stage, we can't be sure of this until a model is found that works. 
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5. OUTLOOK 

The first two years of RHIC running have brought a rich harvest of hadron produc- 
tion data and, as I discussed, abundantly fulfilled our hopes and expectations of finding 
more and stronger evidence for the making of quark-gluon plasma in heavy-ion collisions. 
The observation of large elliptic flow, with its compelling interpretation in terms of fast 
thermalization, and the discovery of jet quenching are two important milestones which go 



significantly beyond what was earlier achieved at the lower SPS energy [ and which 
bring us closer to a generally accepted "proof" of QGP creation. But much more is to 
come: only now, with RHIC finally running at full energy and luminosity (and, hopefully, 
for the full promised time per year) it is possible to address such hallmark measurements 
as thermal dilepton and direct photon emission and heavy quarkonium production, all 
of which play crucial roles in the early diagnostics of the QGP which we are apparently 
mass-producing at RHIC. While trying to solve the HBT puzzle and to quantitatively 
understand jet quenching, we are looking forward to these high-luminosity measurements 
and any surprises they may bring. 

REFERENCES 

1. F. Karsch, Nucl. Phys. A698 (2002) 199c. 

2. J.D. Bjorken, Phys. Rev. D 27 (1983) 140. 

3. R.C. Hwa and K. Kajantie, Phys. Rev. D 32 (1985) 1109. 

4. H. Sorge, Phys. Rev. Lett. 78 (1997) 2309; ibid. 82 (1999) 2048. 

5. PF. Kolb et al, Nucl. Phys. A696 (2001) 175. 

6. J.Y. Ollitrault, Phys. Rev. D 46 (1992) 229. 

7. K. Adcox et al. (PHENIX Collab.), Phys. Rev. Lett. 87 (2001) 052301. 

8. D. Molnar and M. Gyulassy, Nucl. Phys. A697 (2002) 495. 

9. U. Heinz and PF. Kolb, Nucl. Phys. A702 (2002) 269. 

10. PF. Kolb, J. Sollfrank, and U. Heinz, Phys. Rev. C 62 (2000) 054909. 

11. K. S. Lee, U. Heinz, and E. Schnedermann, Z. Phys. C 48 (1990) 525; E. Schneder- 
mann, J. Sollfrank, and U. Heinz, Phys. Rev. C 48 (1993) 2462. 

12. M. van Leeuwen et al. (NA49 Collab.), |nucl-ex/ 0208011 

13. F. Antinori et al. (WA97 Collab.), Eur. Phys. J. C 14 (2000) 633. 

14. H. van Hecke, H. Sorge, and N. Xu, Phys. Rev. Lett. 81 (1998) 5764. 

15. N. Xu and M. Kaneta, Nucl. Phys. A698 (2002) 306c. 

16. J.M. Burward-Hoy (PHENIX Collab.), |nucl-ex/0206UT^ and [nucl-ex/02lOT?01 . 



17. T. Peitzmann, Eur. Phys. J. C, in press ||nucl-th / 02070 1 2 



18. See talks by G. Van Buren ||nucl-ex/ 0211112111 and C. Suire |frucl-ex/021H7T7fl in the 



Quark Matter 2002 proceedings, Nucl. Phys. A, in press. 

19. P. Huovinen et al., Phys. Lett. B 503 (2001) 58. 

20. U. Heinz and PF. Kolb, |hep-ph/020406T| 

21. U. Heinz, Nucl. Phys. A685 (2001) 414. 

22. K.H. Ackermann et al. (STAR Collab.), Phys. Rev. Lett. 86 (2001) 402; C. Adler et 
al. (STAR Collab.), Phys. Rev. Lett. 87 (2001) 182301; ibid. 89 (2002) 132301; Phys. 
Rev. C 66 (2002) 034904; and |nucl-ex/0206U0^| . 



23. K. Adcox et al. (PHENIX Collab.), Phys. Rev. Lett. 89 (2002) 212301. 



10 



24. S. Esumi et al. (PHENIX Collab.), |nucl-ex/ 02 10012 



25. B.B. Back et al. (PHOBOS Collab.), Phys. Rev. Lett. 89 (2002) 222301. 

26. PF. Kolb et al., Phys. Lett. B 500 (2001) 232. 

27. D. Teaney, J. Lauret, E.V. Shuryak, Nucl. Phys. A 698 (2002) 479; |nucl-th/0110037 . 

28. S. A. Bass and A. Dumitru, Phys. Rev. C 61 (2000) 064909. 

29. T. Hirano and K. Tsuda, |nucl-th/ 0205041 - 

30. Z.W. Lin and CM. Ko, Phys. Rev. C 65 (2002) 034904. 

31. T.J. Humanic, |nucl-th/0205053| . 

32. R. Baier, A.H. Mueller, D. Schiff, and D.T. Son, Phys. Lett. B 502 (2001) 51. 

33. R. Hagedorn, Nuovo Cim. Suppl. 3 (1965) 147. 

34. C. Slotta, U. Heinz, and J. Sollfrank, AIP Conf. Proc. 340 (1995) 462. 

35. F. Becattini and U. Heinz, Z. Phys. C 76 (1997) 269. 

36. S.A. Bass et al, Phys. Rev. C 60 (1999) 021902. 

37. U. Heinz, Nucl. Phys. A638 (1998) 357c, ibid. A661 (1999) 140c, and J. Phys. G 25 
(1999) 263; R. Stock, Phys. Lett. B 456 (1999) 277. 

38. J. Letessier et al., Phys. Rev. D 51 (1995) 3408. 

39. P. Braun-Munzinger et al., Phys. Lett. B 518 (2001) 41. 

40. P. Braun-Munzinger, I. Heppe and J. Stachel, Phys. Lett. B 465 (1999) 15. 

41. G. Torrieri and J. Rafelski, Phys. Lett. B 509 (2001) 239; J. Rafelski, J. Letessier, and 
G. Torrieri, Phys. Rev. C 64 (2001) 054907 [Erratum ibid. C 65 (2002) 069902]. 

42. P. Fachini et al. (STAR Collab.), |nucl-ex/0211001| . 

43. J.D. Bjorken, FERMILAB-PUB-82-59-THY (unpublished). 

44. R. Baier, |hep-ph/ 0209038 . 

45. R. Baier et al, Phys. Lett. B 345 (1995) 277; Nucl. Phys. B483 (1997) 291; ibid. B484 
(1997) 265. 

46. M. Gyulassy, P. Levai, and I. Vitev, Nucl. Phys. B571 (2000) 197 and Phys. Rev. 
Lett. 85 (2000) 5535. 

47. U.A. Wiedemann, Nucl. Phys. B588 (2000) 303; ibid. A690 (2001) 731. 

48. K. Adcox et al. (PHENIX Collab.), Phys. Rev. Lett. 88 (2002) 022301. 

49. C. Adler et al. (STAR Collab.), Phys. Rev. Lett. 89 (2002) 202301. 

50. S. Mioduszewski et al. (PHENIX Collab.), |nucl-ex/0210"02l| . 

51. I. Vitev and M. Gyulassy, |hep-ph/0208T08] P. Levai et al., Nucl. Phys. A 698 (2002) 
631; X.N. Wang, Phys. Rev. C 61 (2000) 064910. 

52. J. Velkovska. "Hadron production at RHIC - Where is the soft-to-hard transition?", 
talk presented at the 2002 APS-DNP Fall Meeting, October 9-12, 2002, Michigan 
State University, East Lansing, MI. 

53. T. Sakaguchi et al. (PHENIX Collab.), Erucl-ex/0209030|. 



54. C. Adler et al. (STAR Collab.), |nucl-ex/0210033 



55. M. Chiu et al. (PHENIX Collab.), |nucl-ex/ 02 11008 

56. B. Miiller, |nucl-th/0208"038| . 

57. E.V. Shuryak, Phys. Rev. C 66 (2002) 027902. 

58. S.A. Voloshin, |nucl-ex/ 021001! 

59. M. Gyulassy, P. Levai and I. Vitev, Phys. Lett. B 538 (2002) 282. 

60. U. Heinz and M. Jacob, |nucl-th/ 0002042 . 



